Transport and recombination through weakly coupled localized spin pairs in 
semiconductors during coherent spin excitation 
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Semi-analytical predictions for the transients of spin-dependent transport and recombination 
rates through localized states in semiconductors during coherent electron spin excitation are made 
for the case of weakly spin-coupled charge carrier ensembles. The results show that the on-resonant 
Rabi frequency of electrically or optically detected spin-oscillation doubles abruptly as the strength 
of the resonant microwave field 7B1 exceeds the Larmor frequency separation within the pair of 
charge carrier states between which the transport or recombination transition takes place. For the 
case of a Larmor frequency separation of the order of 7B1 and arbitrary excitation frequencies, the 
charge carrier pairs exhibit four different nutation frequencies. From the calculations, a simple set 
of equations for the prediction of these frequencies is derived. 

PACS numbers: 76.30.-v 76.70.Hb 76.90.+d 72.20.-i 
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INTRODUCTION 



Electrically and optically detected magnetic resonance 
experiments (EDMR and ODMR, respectively) are al- 
ternative ways to detect electron spin resonances (ESR) 
in materials with charge carrier transport or recombi- 
nation transitions that are governed by spin-selection 
ruleo 1 ! 2 ! 3 ! 4 ! 5 ! 6 : 7 . The advantage of EDMR and ODMR 
in comparison to the traditional ESR spectroscopy is the 
sensitivity of these methods which is typically 6 to 10 
orders of magnitude higher 3 ' 4 . This has become partic- 
ularly useful for the investigation of paramagnetic cen- 
ters in highly diluted matrices or low dimensional semi- 
conductor thin film devices and interfaces, point de- 
fect o 1 ! 2 ! 8 ! 9 ! 10 and defect clusters. One of the challenges of 
EDMR and ODMR spectroscopy is that the information 
obtained from these experiments is different in compar- 
ison to the ESR data 11 . The reasons for the discrepan- 
cies between ESR and EDMR/ODMR are mainly due to 
the two different measurement approaches which imply 
two different observables: When the density operator p 
represents the spin ensemble to be investigated, the ob- 
servable corresponding to ESR experiments will always 
be spin polarization (P) = Tr(Pp) represented by the 
spin polarization operator P whereas for the indirect de- 
tection through spin-dependent transport or recombina- 
tion, the observables are the permutation-symmetry or - 
antisymmetry operators represented by the singlet 15) (51 
or triplet operators |Tj)(Tj|, respectivelji 1 ^. For many ex- 
perimental EDMR/ODMR studies (the so called contin- 
uous wave experiments) the different description of the 
observables is not relevant, since these experiments are 
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carried out in the incoherent time regime where only a 
line shape analysis of the respective spectra is feasible. 
However, when coherent effects are studied with pulsed 
techniques (pEDMR/pODMR) 13 i 14 i 15 i 16 i 17 the interpre- 
tation of the experiments strongly relies on the proper 
theoretical description of spin interaction during coher- 
ent microwave excitation&ii 

An example for the difference between a pEDMR signal 
and an ESR signal which come from the same spin en- 
semble are weakly exchange and weakly dipolar coupled 
distant pair states in the band gap of an arbitrary semi- 
conductor material with weak spin-orbit coupling as de- 
scribed analytically by Boehme and Lipsii Weak spin- 
spin coupling means the coupling is much weaker than 
the difference Alu = cu a — ujb of the Larmor— frequencies 
Wa,b — 9a - bfJ f ^ Bo of the pair partners a and 6, respec- 
tively. Note that \xb represents Bohr's magneton and Bq 
the magnitude of an external magnetic field to which the 
spin pair is exposed. Examples for such systems could be 
donor-acceptor pairs whose distance is sufficiently large, 
yet not large enough to make donor-acceptor recombina- 
tion impossible 18 , or equivalently, trap-dangling bond re- 
combination in disordered silicon materials such as amor- 
phous or microcrystalline silicon 8 ' 19 . Weak spin-orbit 
coupling is required in order to ensure spin-conservation 
and therefore, a spin-selection rule. It is fulfilled for in- 
stance in all known silicon morphologies but also in many 
organic semiconductor materials 18 . When the two pair 
constituents are manipulated identically with a coherent 
pulse of high field strength^ 3 . B 1 (St^J, ■- ~ / B 1 > Au' 24 , 
whereas 7 is the gyromagnetic ratio), they undergo a si- 
multaneous spin-Rabi oscillation. This can be detected 
by means of pulse length dependence measurements with 
both, pESR and pEDMR. With pESR, the measurement 
could be conducted by integration of the free induction 
decay and would be called a transient nutation experi- 
ment whereas with pEDMR, the rate relaxation after 
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the coherent excitation would be integrated reflecting the 
pair permutation symmetry within the pairs at the end 
of the exciting pulse 21 . While both, the pESR and the 
pEDMR transients would exhibit oscillating signals, the 
frequency of these oscillations would differ by a factor of 
2: The pESR detected nutation frequency Oesr = 7-Bi 
would simply represent the Rabi frequency of an uncou- 
pled spin s = i, whereas the pEDMR measured oscilla- 
tion would exhibit the frequency at which the identically 
precessing spins of the two pair partners cross the geo- 
metric plane transverse to the field direction of the ex- 
ternally applied magnetic field Bq (the x-y plane) since 
at these moments, the projection of the parallel oriented 
spins in the x-y plane onto the spin eigenstates with sin- 
glet content will be maximized. Since this plane is passed 
twice per nutation period, the oscillation of the transi- 
tion rate is twice as high. Note that this frequency dis- 
crepancy of the oscillations between pESR and pEDMR 
detected transient nutations is changed as the By field 
strength becomes weak: When 7B1 <C Aoj, an ESR exci- 
tation will be possible with only one pair partner at the 
same time. Thus, the maximum singlet content of a pair 
will be achieved when the spin orientations therein point 
into opposite directions parallel to the B field. Hence, 
for weak B\ fields, the nutation frequencies for pESR and 
pEDMR become equal. 

In the following, a semi-analytical study is pre- 
sented which describes spin-dependent electronic tran- 
sition rates (e.g. recombination or hopping transport) 
when the driving forces for Rabi oscillation (under ex- 
perimental conditions these are typically strong coherent 
microwave fields) are between a weak (7B1 <C Alu) and a 
strong excitation regime (7-B1 ^ Aui). Goal of this study 
is to fill the gap between the two analytically derived ex- 
tremal cases of very weak and very strong excitation as 
presented by Boehme and LipsM and to describe a gen- 
eral behavior of the nutation frequency reflected by the 
spin-dependent transition rates for arbitrary B\ and Aw 
and arbitrary excitation frequencies uj. For a straight 
forward interpretation of experimental results 22 , it is of 
particular interest to understand if the change of Rabi 
frequency from Q = f^ESR to fi = 217esr takes place 
continuously or abruptly. 



II. MODEL FOR SPIN-DEPENDENT 
RECOMBINATION 

The basis for the results presented in the following 
are the pair models for spin-dependent recombination 
and transport as described and discussed in detail in 
Refsii^ andSi, respectively. These models are based on 
the Kaplan-Solomon-Mott models under consideration 
of non-negligible triplet transition probabilities and spin 
relaxation. For the calculation of the data presented we 
strictly follow these models under the assumption of neg- 
ligible spin-spin interactions. Note that this constraint 
does not apply to all known experimental systems and 



will always have to be considered when the results pre- 
sented in the following are applied to the interpretation 
of experimental data. 

Quantitatively, the models outlined can be represented 
by an ensemble of spin s = h pairs described by the 
density operator p as derived in RefA We define 
H = Hq + Hi(t) to be the Hamiltonian of an individ- 
ual pair with Ho = —\g a H-BBo<j* — \gb^BBod^ rep- 
resenting the unperturbed Hamiltonian in the presence 
of a constant magnetic field Bo = Bqz and H\{t) — 
— ^9^bBi (<t* + 5+) e~ luJt is the perturbation with a cir- 
cularly polarized microwave of angular frequency u> and 
strength B\. Note the absence of spin-spin coupling in 
Hq in contrast to Eq. 5 of Refiii The dynamics of the 
ensemble of spin pairs can be described by a Liouville 
equation dtp = ^[/3, H]~~ in which, in contrast to Eq. 1 
of Ref, l i, all stochastic terms have been dropped since in- 
coherent processes are considered to be negligible during 
the applied coherent ESR pulse. Only coherent pulses — 
these are pulses that are shorter than the fastest inco- 
herent processes — are considered in the following. The 
ensemble of spins as represented by the density operator 
p = p(t) can be expressed by a unitary, time-dependent 
4x4 matrix. Due to the absence of coupling, the Hamil- 
tonian Ho will be diagonal in the product base with four 
eigenvalues ±-^=p- and ±^A£. wherein loo and Auj repre- 
sent the sum and the difference of the Larmor frequencies 
uj a and Wh, respectively. 

When a solution for p(t) is found, the transients of the 




FIG. 1: Plot of the observable Q as a function of the pulse 
length t and the applied microwave frequency v for a spin pair 
with large Larmor separation 7B1 <C Au). The parameters of 
the simulation were £a = 9.95GHz, = 10.05GHz for < 

Z7T ' ITT — 

t < lps and microwave excitation frequencies of 9.86GHz < 
v = < 10.14GHz as well as a Bj -field of ^ = 10MHz. 

Ait — Z7r 

One can distinguish the two resonant peaks and recognize the 
nutation on the pulse length axis. 
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FIG. 2: (a) Three dimensional color plots of the semi-analytically calculated f2 = FT{Q(r)} as a function of the excitation 
frequency ui scaled in units of 7B1 as the difference between uo and the average of the Larmor frequencies of the two pair 
partners uJ — \(io a + Ub) and the Rabi-nutation frequency fl in units of 7B1. For all four plots (a.i) to (a.iv), tLB\ = 10MHz 
and sr- = 10GHz. From plot (a.i) to (a.iv), the Larmor separation increases: (a.i), 4^ = 1MHz; (a.ii), 4^ = 5MHz; (a.iii), 
^ = 20MHz; (a.iv), 4^ = 40MHz. The two dimensional inset plots display the data of the three dimensional plot at the 
Larmor frequency slices indicated by the dashed lines, (b) The plot of the Rabi-nutation frequencies as obtained from Eqs. [3] 
and ^ for the same parameters as used for the calculated results in (a). A comparison with the frequencies therein shows an 
excellent agreement. 



spin-dependent transition rate 



R(t) = J2nTr[\i}{i\p(t)} 



(1) 



can be calculated from the projection of the permuta- 
tion symmetry operators on the ensemble state. Experi- 
mentally, a real-time observation of R(t) on typical time 
scales in the lower ns-range is difficult to obtain with 
pODMR and often impossible with pEDMR due to the 
insufficient time resolution caused by the long dielectric 
relaxation times of semiconductors in particular at low 
temperatures. Hence, transient nutation experiments are 
typically conducted by means of decay transient mea- 
surements as a function of the applied pulse length 15 . 
The spin-dynamics during the coherent spin excitations 
are obtained from these pulse length dependence mea- 
surements by charge integration Q{t) which reveals the 
permutation symmetry state at the end of the microwave 
pulse as explained in detail in Ref^i. Under the given 
conditions, 



Q(t) oc A(t) 



Pn( T ) 



Pll 



Tr[p s 



Tr[p £ 



(2) 



whereas p^ is the steady state value of the matrix element 
pa of the density matrix p. Hence, it is Q{t) oc A(t) 
which is the observable calculated and displayed in the 
following. Note that while Q(t) represents a number of 
charge carriers for pEDMR experiments, the integration 



of the photoluminescence decay transient in pODMR re- 
veals a number of photons. Nevertheless, in both cases, 
the observable shall be referred to as Q(t) in the follow- 
ing and is always plotted in arbitrary units because of 
this ambiguity. 



III. SIMULATION METHODS AND RESULTS 

When incoherence is negligible, the Hamiltonian H 
can be diagonalized for any given parameter set and 
thus, the Liouville equation in its time integrated 
form p(r) — exp(—j^HT)p s exp(jrHT) can be solved 
by a simple matrix multiplication in which p s = 
\ [\T + )(T + \ + |T_)(T_|] is the same initial state as used 
in Refiii The solution p(r) is referred to as "semi- 
analytic" since this is easy to perform by means of calcu- 
lation of the eigenvalues and eigenvectors of H with any 
given parameter set while attempting to solve it fully an- 
alytically for arbitrary variables leads to unreasonably 
lengthy expressions. 

Figure ^ displays an example for the time domain re- 
sult for Q(t) as a function of the Larmor frequency for a 
spin pair ensemble with a Larmor separation much larger 
than 7S1 which is one of the extremal cases discussed in 
Refii^ (V.A.2.). One can recognize the two well distin- 
guishable peaks at v = 9.95GHz and v = 10.05GHz on 
the Larmor frequency-scale determined by the choice of 
the respective g-values of the spin-pair partners as well 



4 



as the undamped Rabi-nutation on the scale of the pulse 
length r. Outside of the resonances, the signal intensity 
drops while the nutation frequency increases. The latter 
is well known and understood from the expression of the 
Rabi frequency as introduced in Eq. 28 of Ref. 14 . It is 
due to the spin nutation about the residual i?o-ficld ABq 
in the rotating frame which increases with the distance 
of the excitation frequency from the resonance frequency. 
Figure ^ shows that the relevant information contained 
in the calculated transients Q(t) are the frequencies as 
well as the amplitudes of the nutation components. Thus, 
from the solutions of Q(t) the absolute Fourier transform 
(FT) was calculated in order to analyze the various nu- 
tation frequencies contained therein. In contrast to the 
short 1/is transient shown in Fig.^ the time scale for the 
time domain simulation and, therefore, the Fourier inte- 
gration was chosen to be bfis long. For the parameters 
used throughout this study this corresponds to at least 
50 Rabi oscillation periods. 

Figure EJa) displays Q = FT{Q(r)} obtained by 
semi-analytical calculations for spin pairs with four 
different Larmor separations (4^ = 1MHz, 5Mhz, 
20MHz, 40MHz) for microwave excitation frequencies of 
9.95GHZ < v = £ < 10.05GHZ with a given Si-Held 
of jBj_ = 10MHz whilc the plotg ( a ^ and ( a iv ) ful _ 

fill the extremal cases of small and large Larmor sepa- 
ration, respectively, the plots (a.ii) and (a.iii) describe 
two intermediate cases with jB\ » Aw. Note that the 
scaling of the color code was normalized to the maxi- 
mum for each graph in order to achieve sufficient con- 
trast. In order to be able to compare the four differ- 
ent cases more easily, the Rabi components FT{Q(r)} 
at the resonance frequency of the two spin partners are 
plotted in the respective graphs (red curves) as conven- 
tional two-dimensional plots. Here, the chosen scaling is 
equal for all graphs. The four cases displayed in Fig.|2{a) 
confirm the hyperbolic increase of the Rabi frequency 
fl = \J (7S1) 2 + (lo — Lu a ,b) 2 as the microwave frequency 
is shifted out of resonance. One can deduce from the 
two-dimensional inset plots of Fig.[2Ia.i) and (a.iv), that 
the on-resonance cases show only one frequency compo- 
nent for the two extremal cases, namely Q = jBi for 
large Larmor separation and = 2"/Bi for small Lar- 
mor separation. This confirms the analytical results of 
Refii4. Outside of the resonances (u ^ Va,b)> the oscil- 
lation splits into two components for the extremal cases 
and in the general, intermediate cases, there are up to 
four different nutation frequencies. 

In order to illustrate the transition from a single = 
7B1 to a single f2 = 2"fBi frequency component with 
decreasing Larmor separation (or equivalently spoken, 
for an increasing microwave field i?i), a plot of Q = 
FT{Q(t)} versus the ratio on a logarithmic scale 

around Ag- = 1 is shown in Fig. [3] for the one spin on 
resonant cases lo = oj a j, (a) and the average frequency 
case lo = lo' = \(oo a + Ub) (b). Again, for the extremal 
cases of Aj$f- <C 1 and AgL ^> 1, the plots confirm the 
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FIG. 3: Color plot of the Rabi-frequency components Q, — 
FT{Q(r)} as function of log ^^q-^) for an excitation fre- 
quency (a) lo = uj a ,b on resonance with one of the pair partners 
and (b) lo = lo' on resonance with the average lo' of the pair 
partners Larmor frequencies. 



known results for small and large Larmor separation de- 
scribed above. In Fig. CHa), one can see for A^ <C 1 

< 0) that there is only one component with 

A 



7-B1 



7-B1 



(lo, l 

fl = 2^Bi . As Ag^ increases and approaches 1 (log 

approaches 0), this component gradually becomes weaker 
and its frequency increases. This behavior can be under- 
stood by the fact that in Fig. E2a), lo is always equal 
to one of the pair partner resonances. When Alo be- 
comes larger, lo ^ lo' and thus, an increase of the ob- 
served nutation frequency takes place. When A^ « 1 

(log w 0), two new nutation components become 

i) A low frequency component which generally 



7-B1 

visible: ( 

is hard to separate from any given offset in the func- 
tion Q(t) and (ii) one component with f2 = jBi. The 
magnitude of the latter rises from very small values for 
A^ <C 1 and increases asymptotically to a maximum 

value for A^- 3> 1 whereas the two other components 
vanish. A complementary view on these changes is given 
by Fig. EJb) where the Rabi-components for an exci- 
tation frequency lo = lo' are plotted versus For 

■A^ <C 1, plot (a) and (b) agree since they represent 
the same physical situation. When A^ increases and 
Lo a ,b =/= lo' = lo, no low frequency component becomes 
visible. The £1 — jBi component, which also becomes 
visible will increase proportionally to Aw since it is off- 
resonant to the applied microwave frequency lo. 



IV. DISCUSSION 

The simulation of spin-Rabi oscillation as observed 
by pEDMR/pODMR reveals that the doubling of the 
nutation frequency Q = 7B1 to Q = 2^B\ with de- 
creasing Larmor separation {A^) is abrupt which means 
that there is no continuous increase of the oscillation fre- 
quency. Instead, only the magnitudes of the various com- 
ponents change in the intermediate Larmor separation 
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regime about w 1. Here, four nutation components 
become visible which can all become significant at the 
same time. Mathematically, it is obvious that the pairs 
consisting of two s = \ particles with four eigenstates 
will reveal four eigenfrequencies for any given general set 
of parameters. We interpret this behavior for the general 
case as the interplay of the two one-spin systems and 
the one two-spin system by which the spin pairs are de- 
scribed in the two extremal cases discussed above. The 
results of the analytical derivation given in Ref. 14 showed 
that the most significant qualitative change that takes 
place when the Larmor separation changes from infinity 
to zero, is that a transition from a one-spin to a two-spin 
system occurs. For the one-spin system, the spin in reso- 
nance (note that here, only one spin can be in resonance) 
fully determines the oscillation of the pairs permutation 
symmetry whereas for the two-spin system, the permuta- 
tion symmetry is determined by the relative movements 
and spin orientations within the pair ensemble. For large 
Larmor separation, at least one spin will remain fixed in 
its initial state, typically an eigenstate with polarization 
parallel to the i?o-axis. As the excitation frequency is 
changed, the system becomes off-resonant, and two one- 
spin contributions 

tta,b = ^] ^Brf + {uu - uj a , b ) 2 (3) 

appear where the frequency of one increases whereas the 
frequency of the other decreases at the same time. The 
pair still behaves like two individual one-spin systems. 
For small Larmor separation, there are always two spins 
in motion as long as the system is on resonance (there 
is then only one resonance line observed). Hence, since 
the relative spin motion of the two spins within the pair 
will now determine the oscillation of the electronic rate 
transition, the beat frequencies 

o P) „ = n a ± fl b (4) 

of the two one-spin nutation frequencies can be expected. 
Here, Q p and f2„ stand for parallel and antiparallel ori- 
entations, respectively. Figure^b) displays four plots of 
the nutation frequencies obtained with the simple terms 
given in Eqs. and 01 for the same parameter sets used 
for the simulation results displayed in Fig. Ufa). While 
this purely phenomenological description of the nutation 
frequencies can not account for the intensities of the nu- 
tation components — this is the reason why there is no 
grey scale gradient in Fig. [2th) — it nevertheless shows, 
that there is an excellent agreement of Eqs.|3] and 0] with 
the semi-analytically calculated frequency patterns dis- 
played in Fig- Etb). 

While simple quantitative predictions can be made for 
the nutation frequencies, there are no straight forward 
formulas for the prediction of the nutation amplitudes. 



As shown in Fig. Ufa), the intensities of the different 
nutation components can have quite complex microwave 
frequency (u>) dependencies. Well separated lines (case of 
large Larmor separation) exhibit Lorentzian line shapes 
determined by the B\ field due to the power broaden- 
ing as described in Refi 1 ^. For intermediate cases, the 
w-dependence becomes much more complex and a given 
nutation component can exhibit a maximum at its cor- 
responding resonance but also a local minimum at the 
resonance of a non corresponding transition. An exam- 
ple for this behavior is the data of Fig. EJa.iii). At a 
frequency of u — lu' = ±jB\ one can recognize both, a 
maximum of the nutation component at Q = jBi but 
also a minimum of the respective other nutation com- 
ponent at f2 w 2jBi. Qualitatively, this behavior can 
be interpreted by consideration of a four-level system. 
Any of the four levels can undergo first order transitions 
into two different states. For excitations which are out 
of resonance with both transitions, the transition proba- 
bility is small but may not be negligible in the vicinity of 
the resonances. However, when one transition is induced 
resonantly, the transition into the non-resonantly excited 
state is quenched at the same time and the intensity of 
its corresponding nutation component is quenched. Note 
that in spite of this qualitative interpretation, quantita- 
tive predictions have to be made by means of the simu- 
lation methods described above. 



V. SUMMARY AND CONCLUSIONS 

In summary, the response of charge carrier trans- 
port and recombination rates through localized electronic 
states in semiconductors to a coherent manipulation by 
magnetic resonance were simulated as they would be 
expected in pEDMR/pODMR experiments. The tran- 
sient response was calculated with the spin excitation 
present for different excitation fields and frequencies as 
well as different Larmor separations within the pairs. It 
was assumed that exchange and dipolar interaction are 
weak and incoherence due to the electronic transitions 
or spin-relaxation is negligible. The presented data were 
obtained by a semi-analytical simulation method. The 
results show that four qualitatively distinguishable nuta- 
tion processes influence the oscillation of the transition 
rates which reduce to one significant contribution in the 
cases of large and small Larmor separations. Simple em- 
pirical equations for the calculation of these nutation fre- 
quencies have been obtained which match the simulated 
data excellently and a qualitative picture for the inter- 
pretation of the nutation intensities has been discussed. 
The presence of the four nutation processes implies that 
changing the Larmor separation or the applied excitation 
field leads to an abrupt and not continuous change of the 
observed nutation frequencies. 
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